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o, 

^^ I Abstract. We show that with any finite partiaUy ordered set, P, one can associate a 

^N I matrix whose determinant factors nicely. As coroUaries, we obtain a number of results 

Ch ! in the literature about GCD matrices and their relatives. Our main theorem is proved 

3 ' combinatorially using nonintersecting paths in a directed graph. 

oo 

Q . 1. Introduction 

w I Let P denote the positive integers and suppose we are given a subset S = {ai, . . . , a„} 

F^ I of P. The corresponding GCD matrix is (5*) = (sij) where Sij = {ai,aj), the greatest 

c^ ■ common divisor of Oj and aj. We say that 5* is factor closed if given any a^ & S and any 

divisor d\ai then d & S. H. J. S. Smith ^24 proved the following beautiful result about 

the determinant of (S). 

^ ■ Theorem 1.1 (Smith). If S = {ai, . . . , an} is factor closed then 

(1) det(5)=0(ai)---0(a„) 

where is Euler's totient function. D 



in 
in 



O 






Since Smith's pioneering paper, a host of related results have appeared in the literature. 
Q ■ For a survey with references, see the paper of Haukkanen, Wang, and Sillanpaa 1^ . We 

will show that many of these are special cases of a general determinantal identity associ- 
ated with any finite partially ordered set P (see Theorem 12.21 below). This result can be 
^ ' proved by factoring the matrix constructed from P, but we choose to give a combinatorial 

proof based on counting families of nonintersecting paths in acyclic digraphs. This tech- 
nique is due to Lindstrom [20] and independently Karlin ^H]- It was later rediscovered 
rS \ and popularized by Gessel and Viennot [TUl [TT] . We will use the approach in Stembridge's 

c^ ' paper pBj . 

The rest of this article is structured as follows. The next section will be devoted to 
proving our main theorem and giving some preliminary applications. In Section El we 
specialize to the case where P is a meet semilattice. This permits us to invert the sums 
appearing in the general case. Section 0] investigates what happens to the determinant if 
the set S is lower closed (the partially ordered set generalization of being factor closed) 
or meet closed. We end with a section containing comments and an open question. 
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2. The main theorem 

We will first review the digraph machinery which we will need to prove our main 
theorem. Explanations of any undefined terms from graph theory can be found in the 
texts of Harary ^2] or Chartrand and Lesniak [H]. For definitions of concepts about posets 
(partially ordered sets), the reader can consult Stanley's book j23]. 

Let D be a finite acyclic digraph with vertices V and arcs A. Let Rhe a. commutative 
ring with identity and suppose we are given a function (weighting) lo : A ^ R. Then any 
directed path p : wo^i ■ ■ -Vk has an associated weight 

k 

i=l 

We also let 

(^{vo,Vk) = 5]]^(p) 
p 
where the sum is over all directed paths p from vq to Vk- Note that this sum is finite 
because D is finite and acyclic. 

Now suppose we are given two disjoint sets of vertices V = {v[, . . . ,v'^} and V" = 
{v'{, . . . , f"} in V. Consider an n-tuple of directed paths vr = (pi, . . . ,p„) where pi goes 
from v'j^ to ff for alH, 1 < i < n. Then we assign weights to vr and to the pair {V, V") in 
a way analogous to the one used in the preceding paragraph 

n 

uin) = Yl uj{pi) and u}{V', V") = ^ uj{n) 

i=l TV 

where the sum is over all vr where no two of the paths in the n-tuple intersect. 

Finally, given a permutation g in the symmetric group Sn, we define iig = {pi, . . . ,pn) 
to be an n-tuple of directed paths such that pi goes from f^' to u'Vn for all i. So the 
n-tuples considered in the previous paragraph would have g = e where e is the identity 
permutation. If any pair of paths in iig intersect then we call the n-tuple intersecting and 
nonintersecting otherwise. We can now state Stembridge's theorem [2^] about enumerat- 
ing nonintersecting paths. 

Theorem 2.1 (Stembridge). Let D he a finite acyclic digraph with disjoint vertex sets 
V = {v[, . . . , t>^} and V" = {v", . . . , f."} such that any vr^ is intersecting if g ^ e. Let 
(D) = (dij) be the matrix with 

dij =uj{v'i,v'j). 

Then 

det{D)=uj{V',V"). a 

Now let P be a finite poset and consider the incidence algebra I{P, R) of P over R 
which consists of all functions F : P x P ^ R such that F{a, b) = unless a < b. The 



GCD MATRICES, POSETS, AND NONINTERSECTING PATHS 

identity element of I{P, R) is the Kronecker delta function 

S(a b) = I ^ if « = b, 
^ ' '^ [^ otherwise. 

We will also need the zeta function of /(P, R) defined by 

1 if a < 6, 



^0 otherwise. 

The zeta function is invertible and its inverse is called the Mobius function /i of I{P, R). 
In other words, // is the unique function in /(P, R) satisfying 

Yl /^(«>c)= Y^ ^{c,b) =6{a,b). 

a<c<b a<c<b 

Now fix a linear ordering of P which can be used to index the rows and columns of a 
matrix M. Since we will be taking determinants, it will not matter which linear order 
is used and we will merely say that the matrix is indexed by P. We can now state and 
prove our main theorem. 

Theorem 2.2. Let P be a finite poset and let F,G ^ I{P,R). Let {P)fg be the matrix 
indexed by P with entries 

(2) p,, = 5^P(c,a)G'(c,6). 

ceP 

Then 

det{P)FG =11 -^('^) a)G{a, a). 

a£P 

Proof. Construct a digraph D as follows. For the vertices of D, take three copies of the 
elements of P which we denote by P', P", and P'". Now put an arc from a' G P' to 
c"' G P'" if and only if a > c in P. Dually, put an arc from c'" G P"' to b" G P" if and 
only if c < 6. Finally, give weights to the arcs by 

u{a', d") = F{c, a) and uj{c"', b") = G(c, b). 

Consider paths p from a' G P' to b" G P". Clearly all such paths have the form 
p : a', c'", b" where c < a and c < 6 in P. So if we take V = P' and V" = P" then we 
have {D) = {P)fg since 

da'b" = c^(a , b") = Y ^('^' «)<^(c, b) = Pab- 

c<a, c<b 

The last step follows from the fact that a term in pab corresponding to some c not satisfying 
the given inequalities is zero. 

We need to show that D satisfies the hypotheses of Theorem 12.11 By construction, D 
is acyclic. Now take any ra-tuple of paths vr^ which is nonintersecting. We will show that 
this forces g = e. First we claim that the path starting at any a' G P' must then follow the 



ERCAN ALTINISIK, BRUCE E. SAGAN, AND NAIM TUGLU 
b C 




a 
Figure 1. A poset P 

arc to a'" G P'". Suppose this is not true for some a' and let c'" be the next vertex on the 
path. Then we must have a > c in P. Now consider the path of TCg starting at c'. Since it 
cannot intersect the previous path, it must go to some d"' with c> d. Continuing in this 
fashion, we can construct an infinite decreasing chain in P, contradicting the assumption 
that P is finite. So our claim is true. By a dual argument, one can show that the path 
starting at a' must continue from a"' to a" and so f^ = e. Furthermore, we have shown 
that this family vTe is the only nonintersecting path family. So, by the way we have defined 
the weights and Theorem 12.11 

det{D) = uiiTe) = Yl F{a, a)G{a, a). D 

aeP 

We should note that this theorem can also be proved algebraically, either as a corollary 
to a result of D. A. Smith [23, Corollary 2] about the incidence algebra of P, or by directly 
factoring the matrix for P as {P)fg = M^Mq where Mp and Mq are the matrices 
corresponding to the incidence algebra elements F and G, respectively, and t denotes 
transpose. 

As an example of Theorem 12.21 consider the poset P whose Hasse diagram is shown in 
Figure d Then using the linear order a, b, c one obtains 

/ F{a,a)G{a,a) F{a,a)G{a,b) F{a,a)G{a,c) \ 

det F{a,b)G{a,a) F{a,b)G{a,b) + F{b,b)G{b,b) F{a,b)G{a,c) 

\ F{a,c)G{a,a) F{a,c)G{a,b) F{a,c)G{a,c) + F{c,c)G{c,c) J 

= F{a, a)G{a, a)F{b, b)G{b, b)F{c, c)G{c, c) 

As a first application, we note a corollary of Theorem 12.21 which simultaneously gen- 
eralizes a theorem of Apostol [Ij and one of Daniloff j9j. Let f,g be arbitrary functions 
from the poset P to the ring R. Then on substituting F{a,b)f{a) and G{a,b)g{a) for 
F{a,b) and G{a,b), respectively, we immediately have the following result. 

Corollary 2.3. Let P be a finite poset. Let (P) be the matrix indexed by P with entries 

Pab = ^F{c,a)f{c)G{c,b)g{c). 

C&P 
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Then 



det(P) = W F{a, a)f{a)G{a, a)g{a). D 



Qkia) 



aeP 

Now consider the poset defined by using the divisor ordering on P„ = {1,2, .. . ,n}. 
Apostol's theorem is obtained by speciahzing the previous result to the case where P = Pn, 
F{a,b) = ({a,b), G{a,b) = G{h/a) for any function G : Pn ^ R, and g{a) = 1 for all 
a E P. In this case the determinant becomes 

det(P„) = /(l)/(2)---/(n)G(ir. 

With this evaluation in hand, Apostol showed that letting /(a) = a for a E Pn and 
G(a, h) = nipla) where /i is the usual number-theoretic Mobius function gives 

det(c(a, h)) = n\ 

where c{a, b) is Ramanujan's sum. 

To obtain Daniloff's theorem, let P = Pn, g{a) = 1 for all a G Pn, and F{a, b) = 
G{a,b) = Vtkib/a) where 

else. 
Now Corollarv 12.31 becomes Daniloff's result, 

det(P„) = /(l)/(2)---/(n). 

Note that one obtains the same evaluation not just for Vtk, but also for any functions 
F,G : P ^ R such that P(a) = G{a) = 1 for all a e P. 

3. Meet semilattices 

Now suppose that our poset is a meet semilattice L so that every pair of elements 
a,b & L have a greatest lower bound or meet a A b. Note that in this case the sum 
can be restricted to c < a A b. This special case of Theorem 12.21 was discovered by 
Haukkanen, Wang, and Sillanpaa ^1]. They also showed that by further specialization 
one could obtain a theorem of Jager ^Hj concerning a unitary analogue of (0), Smith's 
evaluation 124, of the LCM determinant, as well as a number of other results. 

It would be nice to be able to compute the value of determinants where each entry is a 
single term, as in Smith's original case, rather than a sum. Since Q now consists of a sum 
over all c below a certain element in L, one can use Mobius inversion to accomplish this. 
Thus we can prove the following theorem of Lindstrom [TU] as a corollary to Theorem 12.21 

Theorem 3.1 (Lindstrom). Let L be a finite meet semilattice and let f G I{L,R). Let 
{L)f be the matrix indexed by L with entries 

Lb = f{a A b, a) 
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Then 

(3) det{L)s = \[\y^ii{c,a)f{c,a)\. 

a€L \cGL / 

Proof. Define a function F G I{L, R) by 

c<a 

ii a < b and F{a, b) = otfierwise. By Mobius inversion, tliis is equivalent to 

/(a,6) = 5^F(c,6). 

c<a 

It follows that 

Lb = f{a Ab,a)= J] F(c, a) = J] F{c, a)C(c, b). 

c<aAb cEL 

So applying Theorem 12.21 we obtain 

det(L)j = Yl F{a, a)C(«, «) = H I 5Z ^(^' ")'^*^^' ") ) ' ° 

aG-L aeL \cGL / 

This result can also be found as an exercise in Stanley's book flf)\ Chapter 3, Exercise 
37]. The special case where f{a,b) depends only on a was proved independently by 
Wilf j2n] ■ Smith himself j21| noted that this theorem holds under the further assumption 
that L is a factor closed subset of P. 

To see how Smith's determinant (0) follows from Theorem 13. ![ just let f{a,b) = a for 
all a E S. Then Sij = f{ai A aj, ai) and, since S is factor closed, 

^ H^c, a) f {c, a) = ^iJ,{a/c) c = 0(a). 

c£S c\a 

4. Lower-closed and meet-closed sets 

In work on analogues of Theorem 11.11 there are two conditions that are commonly 
imposed on the set S. If S* C P for some poset P, we say that S is lower closed or a 
lower order ideal if a G S* and b < a in P implies b E S. This corresponds exactly to 5* 
being factor closed if it is a subset of P ordered by division. If, in particular, our poset 
is a meet semilattice L and S is lower closed, then S is also a meet semilattice with the 
same Mobius function as L. So our previous results cover this case without change. 

If S* C L for a meet semilattice one can also talk about 5* being meet closed which 
means that if a, 6 G 5* then aAbES where the meet is taken in L. Again, 5* is also a meet 
semilattice and so Theorem 13.11 still applies and in fact generalizes a result of Bhat j3] 
who considered the case when / is a function of only one argument. 
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However, if 5* C L is meet closed we will show that there is also an analogue of 
which involves the Mobius function of L, rather than that of 5*. This result will generalize 
a theorem of Haukkanen ^3] who, as in Bhat's paper, only looked at functions of a single 
variable. Specializing yet further to L = P„ and the function f{a) = a for all a G Pn, one 
obtains a theorem of Beslin and Ligh pi . 

To state our result, it will be convenient to have some notation. Let £ be a linear 
extension of the partial order on S so that if £ = ai, 02, ..., ctn then a^ < aj implies i < j. 
If (i G -L then we will write d < ai ii d < ai and d ^ aj for any j < i. 

Theorem 4.1. Let L be a finite meet semilattice and let f G I{L, R). Suppose S ^ L is 
meet closed and fix a linear extension £ = ai, 02, . . . , ctn of S . Then 

n / ^ 

det(^)/ = Jjf^ X^/^(c,c?)/(c,ai) 

Proof. We will use unsubscripted variables for elements of L which are not necessarily in 
S. Define F G /(L, R) as in the proof of Theorem 13. II so that we have 

(4) /(ai,aj) = ^F(d,a,). 

d<ai 

Also define F G 1(5, R) by 

F{ai,aj) = ^ ^iJ,{c,d)f{c,aj) = ^F{d,aj) 

d<ai cdL d<ai 

if ttj < aj and F{ai, Oj) = otherwise. If we can show that 

f{ai,aj) = ^ F{ak,aj) 

then the rest of the proof will follow as in the demonstration of Theorem 13.11 So, by the 
definition of F, it suffices to show 

(5) /(ai,aj)=^ ^F{d,aj). 

We will do this by showing that there is a one-to-one correspondence between the terms 
in (0) and those in (JH). 

First note that each d & L occurs at most once in (j3)) and at most once in (jSJ (since 
d<ak for at most one ak E S). li d occurs in (jSJ then d < Ok < Oi and so d occurs in (JD). 
Conversely, if d occurs in (JU then d < ai and so we must have d < Ok for some a^ with 
k < i. But now d < a^ A a^ = ai for some I since S is meet closed. Furthermore, a/ < a^ 
implies I < k, and so I = k since d < a/.. It follows that d < a^ where ak = ai < ai. Thus 
d occurs in (j3)) and we have finished the proof. D 

We should note that Breslin and Ligh [^ have derived a formula for (5*)/ for any subset 
S" C P in terms of an arbitrary lower-closed set containing S. Since this identity expresses 
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{S)f as a product of two (not necessarily square) matrices, one can appy the Cauchy-Binet 
formula to it, as done in Li's paper jTH]! to obtain det(S')j as a sum of determinants. This 
approach was also generalized to meet semi-lattices in an article of Haukkanen |1HJ . 

5. Comments and an open question 

Our results can also yield other information about the matrices {P)fg and {L)f. The 
following theorem is an example of this. In it, we specialize the ring to be the complex 
numbers C for simplicity, although more general rings can be used. 

Theorem 5.1. Let P be a finite poset and let F,G ^ I{P, C). 

(1) The matrix {P)fg 'is invertihle if and only if F{a, a), G{a, a) 7^ for all a E P . 

(2) The matrix {P)fg '^s positive definite if and only if F{a, a)G{a, a) > for all 
aeP. 

Proof. Part (1) follows immediately from Theorem 12.21 and the fact that a matrix over C 
is invertible if and only if its determinant is nonzero. 

For part (2), let the total order used to index the rows and columns of P be a linear 
extension of P. Then each principal submatrix of P is indexed by a lower-closed subset 
S of P. It follows that the submatrix indexed by S is exactly {S)fg- But now we are 
done by Theorem 12.21 again, since a matrix is positive definite if and only if the principal 
subdeterminants are all positive. D 

If L is a finite lattice, then one can compute the determinant of a matrix involving joins 
(least upper bounds) a V 6 of elements a,b E L by working in the dual of L so that joins 
become meets. For general results about join and meet matrices, see the paper of Korkee 
and Haukkanen [T7] . 

However, there are factorizations of such determinants which we have been unable 
to obtain by our methods. As an example, we consider the matrix T„(g) of chromatic 
joins introduced by Birkhoff and Lewis [5j to study the chromatic polynomial of planar 
maps. To define this matrix, let g be a formal parameter. Let H„ denote the lattice 
of partitions of {1,2, ... ,n} ordered by refinement and let NCn denote the lattice of 
noncrossing partitions of the same set. (An excellent survey about noncrossing partitions 
can be found in the article of Simion [221 •) Then Tutte J2H] showed that the matrix of 
chromatic joins could be defined as 

where bk(a) is the number of blocks (subsets) in the partition a. 

Tutte j27j also derived a product formula for the determinant of r„(g) in terms of 
Beraha polynomials. (See also Dahab [Hj.) Letting [-J denote the floor or round down 
function, the nth Beraha polynomial is deflned to be po{(l) = for n = 0, and for n > 1 

ln/2] 






qln/2]- 
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Using the version of Tutte's formula in the paper of Copeland, Schmidt, and Siniion [7] 
gives 



n-l 



detT„(g)=g(''^) JJ 



m=l 



Pm+2{q) 

. QPmiq) 



m + 1 / 2n ^ 
n \n — m—l/ 



The same paper also contains a related determinant-product identity which the authors 
note that they were unable to prove using Lindstrom's Theorem f Theorem 13.11 above). 
although other proofs exist. It would be interesting to find a way to apply the machinery 
of this paper to these identities. 

Acknowledgement. We would like to thank Vic Reiner for helpful discussions and, in 
particular, referring us to the paper of Copeland, Schmidt, and Simion [7j. 
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